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Abstract 

We study new separable orthogonally transitive abelian G2 on £2 models 
with two mutually orthogonal integrable Killing vector fields. For this purpose 
we consider separability of the metric functions in a coordinate system in which 
the velocity vector field of the perfect fluid does not take its canonical form, 
providing thereby solutions which are non-separable in comoving coordinates 
in general. Some interesting general features concerning this class of solutions 
are given. We provide a full classification for these models and present several 
families of explicit solutions with their properties. 

1 Introduction 

The study of spatially inhomogeneous cosmological models has been of great interest 
during past and present decades ||, [|, £|, [j| || . The fact that the Universe is not ex- 
actly spatially homogeneous and the possibility of obtaining more general solutions for 
very-early or late universe models are the main reasons for such studies. The spatially 
inhomogeneous models which have been studied systematically are those space-times 
admitting a maximally 2-dimensional group of local isometries acting on spacelike sur- 
faces (called G2 on S2). Wainwright's classification || for the abelian G2 case is a 
very useful tool in order to deal with such solutions. This classification contains four 
sub-cases depending on the features of the two Killing vector fields that generate the 
group. The most simple sub-case, called B(ii), appears when the Killing vectors are 
both hypersurface-orthogonal (and hence the group acts orthogonally transitively) and 
mutually orthogonal. This class with a perfect fluid source has been already treated 
under some extra conditions, as for example, the imposition of additional homothetic 
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or conformal Killing vector fields (|7|, |J and references therein). The most studied sim- 
plification has been, however, the assumption of separability of the metric functions in 



diagonal and canonical form P, 10, 11 



There have also been general studies on orthogonally transitive G2 cosmologies 
from a qualitative point of view, analyzing the autonomous system of first-order partial 
differential equations coming from the Einstein field equations by using methods from 
the theory of dynamical sistems [|l^, 13, [HJ. The relations between some of the known 



explicit solutions and these theoretical studies are also analyzed in [13] and many 
references therein. 



It must be stressed that the majority of the explicitly known solutions || [10], [LI 
have been obtained by means of the separability of the metric functions in comoving 
coordinates, that is to say, such that the velocity vector of the fluid takes its canonical 
form u oc dx°. This involves two types of restrictions because the separability of 
the metric functions is assumed in a particular well-defined coordinate system. Thus, 
solutions which are separable in other non- comoving coordinates have not been studied 
so far. Notice that the choice of comoving coordinates will destroy in general any 
previously assumed separability of the metric functions. By the way, non-comoving 
coordinates have been already used before (a pioneering paper is [|Tl| ) . 

All this is fully explained in Sections 2 and 4. Section 2 is devoted to presenting 
some general properties of the orthogonally transitive abelian G2 models, in particular 
the so-called "coordinate interchange symmetry" which allows to change the names of 
the non-ignorable coordinates for general G2 diagonal perfect-fluid spacetimes. In fact, 
this property holds also for non-diagonal G 2 models (see ||). Section 3 is devoted to 
showing a broadly general property of spatially inhomogeneous perfect-fluid models, 
that is, that they are in general extendible spacetimes. In other words, solving the 
Einstein field equations for a perfect fluid one obtains regions where the perfect-fluid 
character of the energy-momentum tensor holds, but these regions may not (and in 
general they will not) be a complete spacetime because the algebraic type of the energy- 
momentum tensor cannot hold everywhere. This is a generic property of perfect- 
fluid inhomogeneous spacetimes, and in particular it holds for the diagonal G2 cases. 
This had not been remarked before precisely due to the traditional use of comoving 
coordinates, because the comoving coordinates hold obviously only at the true perfect- 
fluid region. Actually, the use of non-comoving coordinates, as in the present treatment, 
provides natural extensions of the incomplete perfect-fluid regions whenever they are 
extendible. These natural extensions keep the G 2 symmetry of the global spacetime 
but, of course, the energy- momentum tensor cannot keep its perfect-fluid form. The 
possible algebraic types allowed for these extensions are also given in section 3. 

The separation Ansatz in general coordinates is defined in section 4 in a precise 
manner. This leads to a classification of the general separable models depending on the 
values of two natural numbers related to the number of linearly independent functions 
appearing in the metric. The full classification is presented and analyzed, including 
the explicit form of the field equations, which now become simple systems of ordinary 
differential equations. The general kinematical properties of the Weyl tensor, the 
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Petrov type and other relevant quantities are also given. 

Finally, section 5 is devoted to the resolution and the study of the properties of 
part of the cases classified in section 4. Several solutions are explicitly given and many 
of them serve as illustrative examples of the method and their general features. In 
particular, the coordinate interchange symmetry is used to get physically adequate 
solutions. Furthermore, the power of the method is manifested by presenting one of 
the solutions in its canonical comoving coordinates, showing that it would have been 
very difficult indeed to have found it in its comoving form despite the fact that it has 
a perfectly simple expression in non-comoving coordinates. Some concluding remarks 
are given at the end of the paper. 

2 Gi cosmological models 

We are dealing with orthogonally transitive abelian G2 models with two mutually 
orthogonal integrableQ Killing vectors. The matter content is assumed to be a perfect 
fluid, thus there exists a time-like vector field u (velocity vector) such that the energy- 
momentum tensor takes the following form: 



where p and p are the energy density and the pressure of the fluid respectively. The 
existence of two commuting integrable Killing vector fields implies that the velocity 
vector field is orthogonal to them and invariant under the symmetry group. This, in 
turn, implies that u is integrable. From this last result it can be proven a theorem (see 
Wainwright 0) that assures the existence of local coordinates {x a }, already adapted 
to both Killings, in which the metric takes a diagonal form and such that u cx d/dx° 
(canonical form of u), so u oc dx°; that is, there always exist coordinates adapted 
simultaneously to the Killings and the velocity vector so that the metric is diagonal. 

Nevertheless, we can also choose non-adapted coordinates to u. Of course, this will 
be relevant only when some restriction has been imposed on the metric functions in 
a coordinate system. Thus, the assumption of separability allows us to distinguish 
two cases: the particular case when the coordinates that bring the metric functions in 
separate form are chosen to be adapted to the velocity vector, and the general one, 
when these coordinates are not further restricted to be adapted to u. In this last 
generalized case, the velocity vector field of the perfect fluid takes its most general 
form (restricted to be integrable, orthogonal to the Killing vectors and invariant under 
the G2 group) in the coordinates that diagonalize the metric and "separate" the metric 
functions. The line-element for coordinates adapted to both Killing vectors but not 
necessarily to the velocity vector field reads 



T a p = (p + p)u a U(3 + pg a/3 , 



u a u a = — 1, 



(1) 



ds 2 = -F 2 dt 2 + F 2 dx 2 + F 2 (F 2 dy 2 + F 3 ~ 2 dz 2 ), 



(2) 



By "integrable" vector fields we mean hypersurface orthogonal vector fields. 
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where F a = F a (t, x) and the Killings are £ = d/dy, ff = d/dz, while the velocity 1-form 
takes its most general form 

u = u (t,x)0° + u 1 (t,x)O l , (3) 

where {u ) 2 — (wi) 2 = 1 in the orthonormal co-basis 

6° = F dt, O x =F 1 dx, 2 = F 2 1/2 F 3 dy, 0° = F 2 l/2 F 3 1 dz. (4) 

The non-zero components of the Einstein tensor for the metric (0) in this frame are 
Sqoj 5*01) S'lij 5*22 and S 33 , so the Einstein equations for a perfect fluid become 

5*22 = 5*33 , (5) 

Soi = (^oo + 'S'22) (Sn — 'S'22), (6) 

where it must be taken into account that the existence of the perfect fluid is restricted 
to the region defined as follows: 

2S 22 + S 00 - S u ^ 0, sign(2S 22 + S 00 - Sn) = sign(S 00 + S 22 ). (7) 

Note that sign(S o + S22) = sign ( Si 1 — S'22) due to flU) whenever S i 7^ 0. Equations 
(^|) and (|6]) are the necessary conditions for having two spacelike eigenvectors of the 
Einstein tensor with the same eigenvalues (p = S22 = S33) and a third eigenvector 
with the same eigenvalue (p). Conditions (0) assure both the existence of this third 
eigenvector as well as its spacelike character. (We refer to section |3] for a deeper 
discussion on this topic.) 

The perfect-fluid quantities defined in region (|7|) are then 

P — S22, P = S22 + Soo — Sn, 

/ n2 'Soo + S22 / n 2 Sn — S 2 2 ,„s 

M = , («i) = , (8) 

p+p p+p 

where the signs for uq and U\ must be chosen such that the relation 

S01 = (p + p)« «i> ( 9 ) 

which comes from (|l|) and Einstein's field equations in units with 8irG — c — 1, holds. 
It still remains, of course, the freedom on the whole sign for u, which is usually chosen 
such that u° is positive (i.e. u is negative). This means that both u and d/dx° are 
future-directed (say). 

Let us introduce now an interesting property of coordinate interchange symmetry 
of the diagonal G2 geometries which will be very useful for the purposes of classifica- 
tion and the obtaining of solutions in the next section. In fact, this property can be 
generalized to non-diagonal G 2 B(i) models (see ||). 
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Coordinate interchange symmetry 

Given any explicit function f(t,x), let us define 

f(t,x) = f(x,t). 

Upon this definition, we can construct a new line element (denoted by ds 2 ) which will 
consist of interchanging the t and x variables in the original metric functions, that is 

ds 2 = g^pdx a dx 13 . 

We can compute the Einstein tensor S a p of the new metric g^p. It is then easy to find 
the relation between S a p and S a/ 3, which is simply: 

5*22 = ^^22, 5*33 — ^>33, 

o _ Q o o o o 

joi — JOl, >->oo — Jll — 

The first consequence of these relations is that the Einstein equations (||) and (|^) are 
invariant under this coordinate interchange in the sense that the Einstein equations 
for the new line-element are the old equations (||) and (f)|) with t and x interchanged 
in the metric functions and, of course, in the derivative operators (changing dots and 
primes). In other words and explicitly 



5*22 — S33 — S22 — S33 — 



(SqO + S22)(S\l — S22) — S 01 •<=>- (Sqo + 5 , 22)(5'ii — S22) — 5*01 . 

In short, this means that if ds 2 is solution of equations fl5|) and (^), so is ds 2 . However, 
notice that the region of existence of the perfect fluid for ds 2 , given by 

25 22 + S 00 - §u ^ 0, sign(2S i 22 + S 00 - S u ) = sign(5 00 + S 22 ), (10) 

does not coincide in general neither with the previous region (0) nor with the "tilded" 
region (0) (with t and x interchanged). This region must be found separately in each 
particular case. In fact, the only thing that can be said in general is that the region 
where 2S22 + S 00 — Su = is the region where 2S22 + Soo — Su = with t and x 
interchanged because 



25*22 + Sqo — Su — — (25*22 + 5qo — 5 



11/ 



With regard to the fluid quantities, the coordinate interchange transformation induces 
the following changes 

(U ) 2 = - (=«-(tf)2»), (u 1 ) 2 = - S ° + S22 (="- W"), (11) 

p + p V ' p + p V ' 

P=~P, P= -P- 

Note that in (0) there is no contradiction in the signs of the quoted expressions because 
they stand for their expressions (without the minus sign) in (^), which are positive in 
the region (0), but change to be negative when they get transformed and defined in 
the new region (10). 
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3 A broadly general property of the G2 solutions 2 



We will assume that both Killing vector fields are globally defined (well defined over the 
whole space-time manifold). If this did not happen, we could take the open subset of 
the original manifold admitting two spacelike isometries as the manifold itself. Notice 
that (0) involves t and x in general, so that (^) restricts the allowed values of the 
coordinates (apart from very particular cases which involve only the parametres of 
the solution). Then, formulae (0) show that there appear different regions for G 2 on 
52 spacetimes depending on the algebraic type of the Einstein tensor at their points. 
Actually, this is a general feature as will be shown in what follows. 

Equations (g) and (§) restrict the Einstein tensor, and hence the energy-momentum 
tensor via the Einstein equations, to have four possible algebraic types: {1,(111)}, 
{(1,11)1}, {(1,111)}, and {(2,11)} in Segre's notation (see JL5] and references therein). 
Only the first three types admit a timelike eigenvector (which is unique only in the first 
case), while in the fourth type there exists a null eigendirection. The explicit canonical 
form of S a/ 3 for each Segre type together with the regions where the respective types 
hold are given below: 

a) S a p = pv°Vp + p (vlpp + v%Up + v^Vp), at the region defined by 

{25*22 + 5 00 - Sn ^ 0, sign(25 22 + 5 00 - S u ) = sign(5 00 + 5 22 )}, 

b) S a p = —a w^Wp + bw^Wp + a (w^Wp + w^Wp), at the region given by 

{25 22 + 5 00 - 5n ^ 0, sign(25 22 + 5 00 - 5 U ) = -sign(5 00 + 5 22 )}, 

ci) S a p oc g a p, at the region with {25 22 + 5 00 - 5n = 0, S 01 = 0}, 
cii) S a p = ck a kp — dg a p, {k a k a = 0, c ^ 0), at the region 

{25 22 + 5 00 - S n = 0, 5 i ^ 0}, 



where both {v a } and {w a } are orthonormal cobases. The case a) corresponds to the 
perfect fluid region which we will call the ^4-region from now on. The zone where the 
case b) holds will be called the £>-region. Here, there is a timelike eigenvector of S a p 
with eigenvalue equal to p = 5 2 2 = 5 33 (in A, all three eigenvectors with eigenvalue 
p are spacelike). The region defined by 25 2 2 + S 00 — Sn = is the border T (cases 
ci) and cii)), that divides the space-time manifold in the A and B regions in the sense 
that for every continuous curve containing points in A and B, there always exists at 
least one point of the curve in T (this is why we call it a border; see also fL6|). The 
assumption of analyticity of the Einstein tensor on the whole manifold assures that T 

2 We have entitled this section referring to the G2 inhomogeneous models which are being treated 
in the present work, but this property or similar ones apply to more general situations. 
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is either the entire spacetime or its interior is empty in the manifold topology, while 
the mere assumption of smoothness implies that the A and £>-regions are open sets. 

The behaviour of the fluid velocity vector (defined in A) when approaching the 
border T varies depending on whether case ci) or cii) holds at T (this can be easily 
seen from (||) and (||), or also from (|)): if Sq\ 7^ at J 7 , the components of u in the 
orthonormal co-basis must diverge when approaching JF, but this does not necessarily 
happen when S i = at points of T. Then, in a general situation there will be parts 
of the border T where u diverges (see [1? for a general discussion). 



Therefore, in general, a space-time corresponding to a solution of equations (|5|) and 
(H) will be divided into three regions (not necessarily connected) depending on the 
Segre type of the energy-momentum tensor at their points. This could be somewhat 
expected due to the quite general conditions we are using. Our main interest, of course, 
will be the perfect-fluid region A, where we can construct a coordinate system in which 
u takes its canonical (comoving) form while keeping the diagonal form of the metric (see 
the previous section 0). These coordinates are only defined in this region in general, 
because the change of coordinates is not valid where u diverges, that is, on JF, so the 
search of perfect-fluid solutions in comoving coordinates may lead to solutions with 
coordinate singularities that would correspond to an T border. 

Therefore, it seems a natural feature of the G2 perfect-fluid solutions to be ex- 
tendible across .F, and the extensions cannot keep the perfect-fluid character over the 
whole spacetime. Of course, the above problem has a solution if we do not use comoving 
coordinates, which is our purpose in this paper. By using non-comoving coordinates 
and solving Einstein's equations (^) and (|]) we obtain, in one single stroke, both re- 
gions A and B and the border T . Thus, by using this method of obtaining perfect-fluid 
solutions we also get extensions of the A regions, that is, of the perfect-fluid regions 
which are extendible across the border T where the (comoving) coordinate singularity 
appears. 



4 Separability in general: A classification for the 
general models 

Due to the diagonal form of the metric (0), the notion of separability will be applied 
to the metric functions in the obvious way: 

F a (t,x) = T a (t)X a (x). 

Once this is assumed, two different cases appear: (i) t,x adapted to u (comoving 
coordinates) and (ii) t,x non-adapted to u (non-comoving coordinates). 

All the possible solutions in case (i) have been already identified and studied in 
|9|, [K], pf . Ref. |J finds the general solution under the extra assumption F 3 = F 3 (t) 
(see (0)) which implies that the three-slices orthogonal to the fluid congruence are 
conformaly flat. In ref. [0, the remaining solutions are identified unless in the special 
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case of p = p. Finally, [O] , provides the p = p solutions which complete the comoving 
case. In this paper, we study the second (non-comoving) case. 
With the help of separability we get the following line element 

ds 2 = e 2/l (-dt 2 + dx 2 ) + e f2 (e 2f3 dy 2 + e~ 2f:i dz 2 ), (12) 

where f a (t,x) = T a {t) + X a (x) (a = 1,2,3), while u takes its most general (non- 
canonical) form (Q). The explicit components of the Einstein tensor in the orthogonal 
frame (|) are given in Appendix [A]. Equations © and © read then 

T 3 + T 2 T 3 = A7 + Xp^ = K, (13) 
f x X' 2 + T 2 (x[ - lx' 2 ) - 2f 3 X'X = (M 7 (t) + N 8 (x))(M 8 (t) + N 7 (x)), (14) 



2 

where the dots and primes denote derivatives with respect to t and x respectively, 
K is a separation constant, and the M(t)'s and iV(x)'s functions stand for definite 
combinations of the first and second derivatives of the metric functions which are 
explicitly defined in Appendix Equation ( |T^D can be re-written in the form (see 
Appendix [A]) 

8 

Mi(t)Ni(x) + M 7 (t)M 8 {t) + N 7 (x)N 8 (x) = 0, (15) 

i=i 

and differentiating this last expression with respect to t and x we get the following 
equation: 

E^v7 = o- (16) 
i=i 

This implies that if we define n and q to be the number of linearly independent functions 
among the sets {M^} and {N-} respectively, we have q < 8 — n, that is, there are, at 
most, 8 — n linearly independent functions among the {N-}. Thus, at first sight, n 
might take all values from to 8, but this will not be necessary eventually because of 
the property of coordinate interchange symmetry we have shown previously in section 
0: thanks to the t «-> x "symmetry", we have a way to relate solutions of equations 



(|i~3|) and (|T4]) interchanging their n and q numbers. Thus, if a solution has n = 3 and 
q = 5, for instance, then by using the coordinate interchange we get another solution 
with q = 3 and n = 5, and so on. Therefore we must only treat the Einstein equations 
for the values of n ranging from to 4. The rest of cases can be studied by means of 
the coordinate interchange. 

Notice however that, in fact, we have only three original functions of t, {T a }. There- 
fore, we can define another integer m as the number of linearly independent functions 
among the {T a } (obviously m runs from 1 to 3 because m = avoids any dependence 
on t of the metric). Of course, m will be related with the previously defined n. The 
combination of both n and m gives us a way of obtaining solutions and a classification 
for them. 
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To proceed with this classification we will start with any value of m, and then this 
m is related on each case (using some results given in the Appendix |B|) with the number 
of linearly independent functions among the set {l,Mj}, which is exactly n+1 due 
to Lemmas 1 and 2 of Appendix ||. Then, each of the three cases m = 1, 2, 3 will be 
divided into the possible values that n can take for each case. For a given pair {m, n} 
some subdivisions may appear depending on the posible relations between the set of 
functions {l,Mj}. These different possibilities will be given below. Nevertheless, at 
this stage, we have not used yet the equations (|T3|) and ([LI]) explicitly, from where 
new relations arise. These new relations together with the previous ones provide the 
systems of differential equations for the t-functions and for the x-functions, from where 
new restrictions may appear after their compatibilization (if necessary). Therefore, 
the aim of the next subsections is to give the steps in dividing the cases for a given 
m (subsections [O], |4[^ and fO|) . Then, in subsetion |4~4] , equations fll3|) and fll4[) are 
treated in order to present the complete set of equations for the t- and x— functions. 
All throughout the next subsections, the latin lower-case characters will stand for 
constants. 



4.1 Case m = 1 

We have for this case T a (t) = c a T(t) and we need to impose T(t) ^ and at least a non- 
vanishing c a to assure one linearly independent function among the {T a }. Furthermore, 
T will be non-zero because otherwise the metric would be static. It is very easy now to 
establish the possible values of n because, for L running from 1 to 6, we have Ml oc T 2 , 
and the only other possible independent function contained in {Mi} is T, appearing in 
M 7 and Mg. Therefore there are at most 3 linearly independent functions involved in 
{l,Mj} and, consequently, n can take the values 0,1, and 2, dividing the case m = 1 
into three subclasses. We give now the characterization of these subclasses for a given 
n that will be used for the full analisys of the case m = 1 in Section |5|. 

I. n — : this means Mj = for (i = 1 . . .8), which is equivalent to T = 1 (the 

constants are absorved by the c a coefficients). 

II. n = 1 : to avoid the previous case the two linearly independent functions must be 

chosen as {1, T 2 } and we must also have M 7 = c n T 2 + b 7 and M 8 = c 81 T 2 + b 8 . 
Nevertheless, b 7 and b$ can be absorved by M 7 and M 8 (see Appendix |A|). This 
case, taking into account the t-equation in (^) and the fact that there is at 
least a non- vanishing c a , implies an equation for T(t) with the following form: 
f = af 2 + b. 

III. n = 2 : now, there is a relation of the form a + bT 2 + cM 7 + dM$ = and there 

appear two different possibilities depending on wether c ^ or c = 0. 

(i) c 0. Then we can take {l,T 2 ,Mg} as linearly independent functions, and 
M 7 = c 71 f 2 + c 72 M 8 . 
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(ii) c = 0(=>- d 7^ 0) Then {l,T 2 ,Mf} can be chosen as linearly independent 
functions, and M 8 = c$iT 2 . 

Again, the possible additional constants have been absorved in M 7 and M 8 as in 
the previous case. 



4.2 Case m = 2 

In this case we have T a (t) = c a T(t) + d a K(t) where {T(t),K(t)} are two linearly 
independent functions. We also need to impose, of course, that the matrix composed 
by c a 's and d a 's has rank two. This case can also be treated dealing directly with the 
{T a } functions and dividing this class depending of which pair is taken to be linearly 
independent, but we have preferred to introduce the functions T(t) and K(t) for the 
sake of compactness and brevity. 

In fact, 1, T(t), K(t) can be assumed to be three linearly independent functions, as 
otherwise we could reduce this class to the case m = 1. For, suppose, on the contrary, 
that there existed a linear relation aT{t) + bK(t) + c = 0. Then, if a = we would 
have that K(t) is a constant and therefore, it could be set equal to zero because the 
terms e daK in the metric can be absorved into the coordinates. If a ^ we would 
have T{t) = —{b/a)K{t) — c/a, thus T a (t) = (—c a b/a + d a )K{t) — c a c/a, and again, 
redefining the constants and absorving the constant terms into the coordinates, we 
could set T(t) = 0, in contradiction. From Lemma 1 it follows then that T and K 
are two linearly independent functions, and therefore, Lemma 3 of the Appendix [E| 
implies that the set {T 2 , K 2 ,TK} consists of 3 linearly independent functions. This, 
in turn, means that among {Ml} (L = 1 . . . 6) there are exactly 3 linearly independent 
functions (see Appendix |FJ) . 

In M7 and M 8 there appear two other functions (T and K) that can be linearly 
independent from the rest. Therefore, in the set {l,Mj} we have a minimum of 3 
linearly independent functions and a maximum of 6 so that n can take the values 2,3,4, 
and 5, but this last case (n = 5) does not need to be treated thanks to the coordinate 
interchange symmetry t <-> x. In summary, the case m = 2 is divided into three 
subclasses n = 2, n = 3 and n = 4. 



4.3 Case m = 3 

In this case the three functions {T a } are linearly independent. Nevertheless, we will use 
here three generic independent functions in order to give a more compact presentation 
of the cases. Thus, let T(t), K(t), Q(t) be three linearly independent functions such that 
T a (t) = c a T(t) + d a K(t) + e a Q(t), where the determinant of the 3x3 matrix composed 
by the constants must be non-zero. 

Analogously to what is explained in the second paragraph of subsection [4.2| , the set 
{l,T(t), K(t),Q(t)} consists of four linearly independent functions, as otherwise we 
could reduce this class to the case m — 2. From Lemma 1 it follows then that T, K, 
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and Q are three linearly independent functions, and Theorem 1 implies then that there 
are at least 5 linearly independent functions among the set {T 2 , K 2 , Q 2 , TK, TQ, KQ}, 
that is, among the Mi functions, and hence, among {1, Mi}. Therefore, n can take the 
values from 4 to 8, so that only the case n = 4 has to be treated. 

The three cases with m = 2 and the single case with m = 3 will not be solved 
explicitly in this paper. Nevertheless, we present the whole set of equations and the 
kinematical quantities in general in the next two subsections. 



4.4 The complete set of equations 

Now, we proceed with the obtention of the complete set of equations coming from 
equations ( |i~3]) and flU] ) once we have chosen the pair {m, n}. As was explained above, 
these equations together with the systems given in the previous subsections will form 
the full set of equations for the t- and x- functions. In fact, equations flT3| ) need no 
further treatment and will not be repeated in this subsection. Thus, we focus on 
equation ( [14] ) written in its form flT5p. 



We choose {1, m^(i)} (A = l...n) to be the n+1 given linearly independent 
functions such that 

n 

Mi(t) = E c iA m A (t) + h {i = 1 . . . 8), (17) 

A=l 

where the constants qa form a matrix of rang n.f] Therefore, we take the functions 
rriAit) and the constants bi to be the fundamental objects with regard to the Mj func- 
tions. Now, from (|16|), using (|i"7|) and the fact that {rriA_{t)} are lineraly independent 



(Lemma 1), there appear n linearly independent relations between N-. This implies 
that there are at most 8 — n linearly independent functions in {N[}. Thus, we choose 
{1,ub(x)} (B = 1 ... (8 — n)) to be the 9 — n = (8 — n) + 1 linearly independent 
functions such that 



8-n 



Ni(x) = V d iB n B (x) + di (i = 1. . 



B=l 



where the constants diB have no restriction a priori. 

Now, derivating equation ([15]) with respect to x and using ([[7]) we get 



E E ^aN[ m A (t) + E hN[ + (N' 7 Ni)' = 0. 



A=l \i=l / i=l 



Using the fact that {l,rriA{t)} are linearly independent functions, from this last ex- 
pression it follows that Y% c iaN- = and J2i hN- + (N 7 N 8 ) = 0. Integrating this pair 

3 Note that we could set 67 = 6 8 = without loss of generality (see Appendix 0) , but this will not 
be always used. 
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of expressions we get 

8 

Y,c iA N l = C A (A=l...n), (19) 

i=l 

8 

Y,kN t + N 7 N 8 = B, (20) 
i=\ 

where C_a and B are constants. Differentiating equation ( |T4| ) with respect to t and 
using ( |18D we obtain the corresponding expressions for the Mi functions. The equation 
analogous to ([TP]) does not give relevant information (only relations between constants 
that will not be used at the end), while the analogous to fl2~0|) reads 

8 

Y^nM, ■ M;,\L A. (21) 
i=i 

We take now the original Eq.flTo]) and substitute the functions Mi from ( |T7| ) and the 
terms M 7 M$ and N 7 N$ isolated from ([H]) and (|20|) respectively. Using again the fact 
that {l,m A (t)} are lineraly independent and taking ( |19|) into account, the following 
relations arise 

8 8 

Y a ibi = A + B, Y a i C iA = CA- 

i=l i=l 

Putting these last expressions into (^T|) we finally get 

n 

Y C A m A + M 7 M 8 + B = 0. (22) 

A=l 

Thus, equation (|T5D splits into the equivalent set of ordinary differential equations (|H|), 
(p0|), and fl2~2"|) . Therefore, the complete set of equations consists of (|13|),(|T9| ) ,(|20| ) , and 
([22]) together with the relations between the t-functions arising from the classification 
of the previous subsections. 



4.5 Kinematical quantities and the Weyl tensor 

For the sake of compactness in the expressions of the kinematical properties of u we 
define first the following objects 

w = e 2 ^(S 00 + S 22 ) = M 7 (t)+N 8 (x), 
Wl ee e 2h (S u -S 22 ) = M 8 (t) + N 7 (x), 

so the Einstein equation (|]) reads simply WqWj^ = S 2 , and the perfect-fluid quantities 
take then the following form 



p + p = e - 2h {w - Wl ), p-p= e~ 2h {f 2 + T 2 - X'i - X' 2 



2 h 



W \2 _ W l 



( u o = > \ u i) ~ 
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The .A-region defined in section |3] is given now by the conditions wo — w\ ^ 0, sign(w — 
w-j) = sign(u> ), which can be combined with the extra requirement p + p > to give 



w > 0, 



w — w 1 > 0, 



(23) 



which become the necessary and sufficient conditions to have a perfect-fluid source 
satisfying the energy condition p + p > 0. From now on, the domain defined by these 
conditions (|23|) will be referred to as the A^-region. The relative signs of u and u x 
are determined then by S = (w — w 1 )u u 1 . 

The fact that u is integrable implies that its vorticity vanishes, that is, u> af 3 = 0. 
The expansion reads as follows, 



= —. V2 { K - ™i) [S (X[ +X' 2 )-w (T, + T 2 

u {w - w x y v ' 

+w S' - w' E + -(w w 1 - WqWi) I . 
The non-zero components of the acceleration computed in the co-basis given in (|J) are 

Mi 



1 . 
u 



a — ^ a l5 



-/i 



(w - wj- 



+w S - w S + -{w' Q w x - w w[] 
while the shear tensor has the following non-zero components 

2 



00 



CT 



Ml 



01 



-a 



"o 



Hi 



11; 



Lin 



a 



11 



(w - Wj) V 



(E^- Wo T 2 ) + -K) 2 ^, 



(7 



22 



W (W - tUiJ 



W T 3 - EX3 



3 

1 1 



2 to,) 2 * 11 ' 



a 



33 



M (W - Wj) 

and the shear scalar 2a 2 = o" Q/ 3<t q/3 is then 



2K) 



a' 



tic 



-9 - e~ h . 

3 [w - v 



I T>X 2 — w T 2 



u (w - w x ) 



l 2 



w T 3 - SX^ 
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The non-vanishing scalars of the Weyl tensor computed in the null tetrad k = 2 1 / 2 (0°- 
1 ), / = 2- 1 / 2 (0° + 1 ), m = 2- 1 / 2 (0 2 + i0 3 ), where (a) are given in (|), are 



+ ^4 = 2e~ 2/l T t T 3 + X[X' 3 - A 



^2 



^ 4 = e 
1 



-2/i 



2 fcx't + X(T 3 ) - T 3 X^ - X' 3 f 2 



2./i 
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4 (T 3 2 - X 3 2 ) + 2 (X'l - Ti) + T 2 - X 2 ' 



(24) 
(25) 

(26) 



where we have only used equations (|T3|) from which we have isolated both T 3 and X 3 
(which appear in \l/o and ^4). From fl24j) -(p6"l) it follows that the solutions will be in 
general (and at generic points) of Petrov type I. 



5 Full analysis of the case m = 1 and explicit solu- 
tions 



In this section we treat the case m = 1 given in subsection LI so that T a {t) = c a T(t). 



We give the equations for each of the subcases taking into account the results of subsec- 



tion |4.4j and then some particular examples will be solved. Unless otherwise is stated, 
throughout this section the generic term "solutions" will stand for maximally G2 not 
included in the previous works on separable comoving coordinates 0, [TI| |IIfl , that is, 
they will non-separable in comoving coordinates "a priori" . Furthermore, a relationship 
between n and the type of equation of state appears. 

5.1 n = 

As was shown in subsection O, we have now T a (t) = c a , so that from fl4"2"|) and (H3) of 



the Appendix |A] and ( |17D we get the values for b; L = {c 2 , c|, c 2 , Cic 2 , Cic 3 , c 2 c 3 , C\c 2 - 
2c§, c 2 (ci — l/2c 2 )}. Note that we have not used the fact that one can have b 7 = b 8 = 
so we set a = (3 = in (££|) and ([45]). 

Equations (p~3| ) become simply K = c 2 C3 and 



XZ + X' 2 X^ = c 2 c 3 , (27) 
while Eqs. (0) provide no relations, ( p2|) fixes i? and finally (^) reads 
8 / 1 \ 

2 hNi + N 7 N 8 = -ca ( c a - -c 2 ) ( Cl c 2 - 2c 2 ), (28) 



where Ni are defined in Appendix [A|. Equations (p7|) and (^8|) form a system of two 
first-order ordinary differential equations for the three unknown functions X' a . Thus, 
in general the solutions depend on an arbitrary function, allowing for several further 
Ansatzs in order to find explicit solutions of this system of equations. In particular, 
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this freedom can be used in principle to demand some extra property of the solutions, 
such as particular equations of state. Indeed, it is possible to find solutions with a 
p = 7p equation of state including 7 = 0, that is, dust models. In fact, X a (x) = 
d a X(x) + l a x such that X"(x) 7^ and p = leads to a family of dust solutions that 
can be generalized to give a bigger family of algebraically general dust G 2 models Jl8|] . 
Nevertheless, these solutions have q = 3 (number of linearly independent functions 
among {N-}) so that they will appear, within our scheme, when n = 3 in the m = 2- 
case after a t <-> x change. In this way, in order to avoid any superposition of solutions 
between the different subcases (using the t <-> x, n <-> q "symmetry") and to keep a 
'coherent' classification, we should look only for solutions with the following values of 
q G {0, 5, 6, 7, 8}. The rest of the cases (q = 1, 2, 3, 4) can be left to the study of cases 
(n = 1, 2, 3, 4, q — 0). The cases with q > 5, following the reasoning given in Section^, 
need at least two linearly independent functions among {X a (x)}. The solutions depend 
on the extra assumption which closes the system fl2"7|)-(]2"BD. 

In the remaining possibility q = we have that X' a (x) = k a , where k a are constants. 
These spacetimes always admit a third Killing vector given by 

d d 1 d 
C = ~k\-Q t + c i-q^ + -jV [hc2 - cik 2 + 2(A;iC3 - fc 3 ci)j ^~ 

1 d 

+-z [kic 2 - cik 2 - 2(/cic 3 - fc 3 Ci)J — . 

Equations (p7|) and (^) give two relations on the constants c a and k a which do not 
imply necessarily the appearance of more isometries, thus the resulting solutions belong 
to the tilted Bianchi perfect-fluid models. Nevertheless the study of such solutions are 
out of the aim of the present work, and will be omitted. 



5.2 n=l 

Following the arguments given in the subsection ^72| and looking at (|i~7|) , in this case 
we have: mi(t) = T 2 , bi = and ca = {c\, c?,, c§, cic 2 , C1C3, C2C3, C71, Cgi}, where C71 
and Csi are two new arbitrary constants and 67 and bs have been absorved in M7 and 
Ms- The equations containing T are then (0) for i = 7, 8 and the t-equation in flilf ), 
which read respectively 



(ciC 2 - 2cfj f 2 - Qc 2 + c^jf + a = c 7X t 
ci - \c 2 \ (c 2 t 2 + f ) + p = c 81 T 2 , 



2 

c 3 (c 2 T 2 + T) = K. 



These expressions together with the fact that the c a cannot vanish simultaneously 
imply an equation for T(t) of the form 

f = af 2 + b, 



15 



where a and b are constants. Using this equation in the three previous expressions and 
taking into account that the functions 1 and T 2 are linearly independent, we get the 
six following constraints for the constants: 



c 7 i = cic 2 — 2c s 



-c 2 a, 



Cl 



a 



;C 2 J(c 2 + a), c 3 (c 2 + a)=0 (29) 
\c 2 + Cl ) 6, /? = Qc 2 - ci) 6, K = c 3 b. (30) 



The only remaining equation involving t-functions is (|22|) , that reads now C1T 2 + 
c 7 iC 8 iT 4 + 5 = 0, implying Ci = B = and 



C71C81 = 0. 



(31) 



The equations flTBD, flT9|) and ( PC| ) for the x- functions become respectively 

X^ + X' 3 X' 2 = c 3 b, y t, c a N i = °> N 7 N 8 = 0. 

i=i 

From this last relation and ( |3l"D it follows that only two different cases may appear (if 
we neglect the solutions separable in comoving coordinates): (a) cj\ 7^ {csi = 
0, N 8 = 0, N 7 ^ 0}) and (b) c 7 i = {c 81 ^0,N 8 ^ 0, N 7 = 0}). The final system 
of equations for both subcases can be written in the following compact normal form 
once we have defined < 



x' 2 ^ 



X[ - (1/2)X 2 and ip = X' z and after using (g| 



c 3 b - <pLf 
ip' = (\c 2 + c 1 )b 



+ 2ip)-2<p 



e 2 

+ p (Cl0 + C2V) 

2 c 2 6 



1 2 

— (ci0 + c 2 ip - 2c 3 ip) 



(32) 
(33) 

(34) 



where {e,/C(^ 0)} = {0,c 7 i}, {l,c 8 i} for cases (a) or (b) respectively. 

In case (a) the solution of the constraints ( |2"9"D (now c 8i = 0) give two possibilities: 
(al) c 2 + a 7^ 0, implying then c 2 = 2ci 7^ 0, c 3 = 0, thus /C = c 7i = 2c\{c\ — a) 7^ 0, 
and (a2) a = — c 2 , which gives /C = c 7 i = 2cic 2 + c|/2 — 2c 3 7^ 0. For case (b), the 
constraints imply c 3 = 0, a = 2cic 2 /(c 2 + 2ci) with c 2 + 2ci 7^ (as otherwise it would 
follow the vanishing of the c a ) and K. — c$i — c 2 (ci — c 2 /2)(c 2 + 4ci)/(c 2 + 2ci) 7^ 0. At 
this point it is interesting to give the conditions (53) explicitly, which read 



case (a) c 7 i > 0, 
case (b) c 8 i > 0, 



c 71 T 2 > (ci0 + c 2 ip - 2c 3 <£>Y 
c 2 sl f 2 < (ci0 + c 2 i\) - 2c 3 ipY 



thus, the definition of the .Ag-region (where p + p > holds) gives in general a proper 
determination of a domain in the manifold and also a condition on the constants, 
which is invariant under the change t «-> x. This last feature is important because 
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these "invariant" conditions allow us to prove some statements about the existence of 
solutions under some extra conditions (kind of equation of state, for example) that 
will hold also after changing n <-> q. Using the coordinate interchange, the previ- 
ous conditions become Cn > and c^X' 2 < (c\4> + c 2 ip — 2c 3 </?) 2 (now the functions 
{<j),ip, (f} are redefined with T a {t) replacing X a (x)) in the case (a) and c%\ > and 
c^X' 2 > (ci4> + c 2 ip — 2c 3 <£>) 2 in the case (b). 

In the present case, imposing any further restriction such as equations of state may 
overdetermine the system of equations and give solutions with more isometries. For 
instance, it can be shown that an equation os state of the form p = 7p implies 7 = 1 
(stiff fluid), even for the comoving solutions, as otherwise a third isometry appears 
acting on the original Killing orbits, which became plane (i.e. a plane G 3 on S 2 ). 

The system of equations for the different cases has not been explicitly solved in 
general, but some particular families have been found under some extra restrictions on 
the constants. In fact, when b = in case (a) the system can be completely solved, 
although the solutions will have q < 4. Analogously to what has been explained for the 
case n — 0, the values of q for the representative solutions of this case n — 1 should be 
{0, 1, 5, 6, 7, 8}. We will present now some results concerning solutions with q = 0,1. 



5.2.1 q = 

The solutions with q = (X' a (x) = k a ) correspond to the singular points (<po, ipo, (j>o) of 
the differential system. Let us begin with the case c 3 7^ (=>• a = — C2) (case (a2)), and 
determine then b from fl32|) . In order to keep w l > we must avoid 0o — 0. Therefore, 



equations fl33|) and ([34]) give 



, c 2 + 2ci 

4>o = <Po, 

2c 3 

(c 3 o - c 2 (f ) [2c\ - c 7 i) c 3 o - (2ci + c n ) c 2 y?o = 0, (35) 

where (see above for case (a2)) c n = 2ciC 2 + c 2 ,/2 — 2c\. These solutions do not admit 
a third isometry nor a barotropic equation of state in general. In fact, it can be shown 
that the only barotropic equation of state that these maximally G 2 solutions can have 
is p = p, and those correspond to the cases with c 3 0o — c 2 (fo = 0. Let us focus our 
attention into the rest of the cases. From (^), and taking into account that neither 
2c\ + C71 nor c 2 can vanish in order to keep C71 > 0, we determine (fo, so finally we have 

_ / c 2 + 2cA _ ^3g0 o? ^ 36 ^ 

V 2c 2 / c 2 

where we have defined Q = {2c\ — c^)/{2c\ + C71) (Q < 1). 

At this stage we can compute p — p and realize that it is positive everywhere, so 
the dominant energy condition cannot hold in the perfect-fluid region. Nevertheless, it 
is important to remark here that, within our treatment of the problem, solutions not 
satisfying energy (or other) conditions may still be relevant because these conditions 
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are not invariant under the interchange symmetry t <-> x. Thus, we must always 
check the physical properties both for the explicitly obtained solution as well as for 
its partner solution with t and x interchanged. Otherwise, the full method will not be 
coherent. An illustrative example is given, in fact, by the solutions above, because the 
new solutions obtained by means of the coordinate interchange are well-behaved. In 
order to see this, let us perform the change t <-> x so that the line-element becomes 

X s 2 = e (2^+0o) i +2c 1 X(x) (_ dt 2 + ^ + 

e <P t+c 2 X(x) ^ e 2(<p t+c 3 X(x)) d y2 + e -2(ip () t+c 3 X(x)) dz 2^ ^ ^ 

where ipo and (p$ are given by (|36|) and the function X(x) satisfies the equation 

-12 Q ±2 



X" + c 2 X r - 



C2 



from where four cases arise:[] 



(i) X(x) — — In cosh (^Pq^/Qx^ , x G (-00,00), 

(ii) X(x) = — In sinh (4>o\[Qx^j , x G (0, 00), 

(iii) X(x) = — In cos ($o<J^Qx^ , <fr \f-Qx e (— 7r/2, tt/2), 

1 

(iv) X(x) = — ln(gx), x G (0, 00), 

c 2 



where Q is positive in the first two cases, negative in the third and Q = (cj\ = 2cf) 
in the fourth. In this last case, the solutions admit the timelike homothetic Killing 
vector d/dt (not parallel to the fluid vector), while for the corresponding H «-> x' 
solutions, the homothetic Killing vector becomes d/dx. The solutions do not have a 
barotropic equation of state unless 2ip + <p Q = (implied by c\ = 0), which gives a 
further isometry, so we will demand C\ 7^ in the following. The energy density and 
the pressure read then 



P + P 



-V: 



^0 + \<pQ 



r 2 



X 



12 



P~ P 



-2/i 



2c 71 



2c( + c 7 i 



so in the region Ae {p + p > 0, and remembering C71 > 0), the dominant energy 
conditions are fulfilled. Note that it can be shown that in the case (i) the condition 
p + p > is automatically satisfied, so the region Ae covers the entire spacetime. The 
quantities involved in the rest of the fluid quantities read 



w = on ip + -4>< 



2^) cV 



w-. 



C71X 



r> 



'■-i i ^0 + ^0) —X'. 



4 The case X" = falls into the previous case with n = q = 0. 
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The case when c 3 = needs also 0o — in order to avoid further isometries, 
and e = 1 (case (b)). After straighforward calculations we find, as in the previous 
family, that p — p is positive everywhere thus violating the dominant energy condition. 
However, we can use the coordinate interchange symmetry again and luckily we get 
new solutions which do fulfil the dominant energy conditions in the perfect-fluid region. 
Finally, the metric is thus given by 

Ts 2 = e 2 ^ 1 cos" 9 ' (q 3 x) (-dt 2 + dx 2 ) + cos 9 * (q 3 x) (e 2Lpot dy 2 + e~ 2ipat dz 2 ) , (38) 

where x G (— 7r/g 3 , 7r/g 3 ), we need to impose <p 2 — ip 2 > 0, and we have defined q\ = 
(f 2 o ~ ^o)/ 2( fl ql = (vl - ^o)/(^o + V>o)> and 



<?3 = 2<fl* 



2(y§ + Vo) 



The energy and the pressure are given by 

p + p = q l 2 e-^ co ^{q,x) (2^ 2 tan 2 (g 3 x) - \qtql 
qi \ I 

P-P = -qhle- 2 ^ 1 cos 9 ? (q 3 x) tan 2 (g 3 x) + , 

\ip + % ) 

from where it is evident that, in the perfect-fluid region p + p > 0, the dominant energy 
condition is always satisfied, as claimed previously. There is no barotropic equation of 
state in general, though. For the rest of the fluid quantities we have 

™o = 2^ 2 %tan 2 (g 3 x), w x = ~<fi<fiqj, £ = q^s^o tan(g 3 x). 
Qi A 

These two families of solutions ( j3~^) and (|38"D with q = are algebraically general. 



5.2.2 q = l 



For q = 1 we need X a (x) = d a X(x) such that X"{x) 7^ 0, which is also sufficient, 
because the system of equations (|32"D- (j34D implies necessarily N 7 ,N 8 oc X' 2 . The 
compatibilization of this system gives some relations on the constants: first, we have 
that e = 0, so that the solutions must belong to case (a), in order to avoid comoving 
coordinates. For 6 = 0, the only maximally G2 solutions without p = p equation of 
state are then given, after redefining the constants, by the following line-element 



ds 2 



a .^ t -2(l-M 2 /9) [_ dt 2 + dx 2 + x l-» ^ x v dy 2 + x -u dz : 



(39) 



where q = p 2 + v 2 — 2p — 1, and the corresponding metric after the t <-» x interchange. 
These solutions admit two further conformal Killing vector fields and have no barotropic 
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equation of state in general. These families were found in || (pp.2320), and we refer 
to this reference for further details. 

For b 7^ we need C3 7^ 7^ d 3 as otherwise there would appear a third isometry, 
so that the solutions must belong to case (a2). Then 

d 2 

a = -c 2 , c 2 = -3-C3, c 1 = d 1 = 0, 

where d 2 and c 2 cannot vanish so as to have C71 > 0. The functions T(t) and X(x) 
are given by T = d2C^/d^T 2 + 6 and X" = — d 2 X' 2 + c^b/ds, and the equation of state 
reads p = p + 1d2C^bjd 3 . To assure the dominant energy condition we should first 
impose b = —d 3 fi 2 /(d 2 c 3 ), which turns out to be also sufficient (in the -4.^-region). 
The coordinate interchange t x gives no further solutions for the corresponding Ae- 
regions (which become, in fact, identical to Ae)- Two different families of solutions 
appear. The first is given by 

ds 2 = -dt 2 + dx 2 + cos 1+2u {(ix) cosh 1 " 21 ' {fit) dy 2 + cos 1 ' 2 " '(jjlx) cosh 1+2v {fit)dz 2 , (40) 

where \x and v > are constants. The conditions for the Ae region give v < 1/2 and 

tanh 2 (/zt) > tan 2 (/ix), 

so that we choose t > and 2 G (— 7r/4/i, 7r/4/x). The equation of state and the fluid 
quantities are given by 



p = /i 2 



- — z/ 2 ^ (tanh 2 (/it) — tan 2 (/ix)) + 1 



P = P~ 2 /A 



tanh (/it) x tan(/ix) 



tanh (/it) — tan 2 (/ix) ytanh (/it) — tan 2 (/ix) 

ai = — /i(-u°-u 1 ) 2 tan _1 (/ix)(tanh 2 (/it) — tan 2 (/ix) — 2), a = — —7-7— raj., 

tanh(/it) 



6 1 = /i (tanh (/it) — tan (/ix) — 1 

'tanh 2 (/it) — tan 2 (/ix) 



/ tan(/ix) \ 2 tan(/ix) 
00 ytanh(/it) J 11 ' " ! tanh(/it) 11 



0"11 



/itanh 2 (/it)[( M °) 2 + (u 1 ) 2 ] / o / \ o / \ 
^2 (^tanh (/it) — tan (/ix) + 2 

3 (tanh 2 (/it) — tan 2 (/ix)^ 



1 



cr 22 = -/iz/y tanh 2 (/it) - tan 2 (/ix) + 2 ^ )2 °" 11 



cr 33 = 2/ivytanh (/it) — tan 2 (/ix) + <r 22 . 
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The non-vanishing components of the Weyl tensor are 

= ^4 = 2zV , ^ 2 = -p 2 - -p, 

thus the Petrov type is I for generic points. 

Notice that this family of perfect-fluid solutions do not present any curvature singu- 
larity. However, this fact may be not relevant because, as explained above, the perfect- 
fluid .A-region is extendible through the spacelike hypersurface tanh 2 (/it) — tan 2 (px) = 
0. A particular extension is in fact given by the line-element (|40| ) itself when taking 
every possible value of t and the range x G (— 7r//x, 7r//x). This particular extension is 
then singular, as can be immediately checked from the above expressions. 

This solution ( f40|) is in fact a good example of the power of using non-comoving 
coordinates. As explained in section 2, every diagonal G 2 solution can be written in 
comoving coordinates. Thus, we can ask ourselves, how does solution ( f40|) look like in 
comoving coordinates? In this case the change of separable coordinates to comoving 
coordinates {£', x'} can be performed explicitly and the line-element (|40| ) in the A- 
region becomes 

ds 2 = 




where f(t',x') = (t' 2 — l) 2 — x' 2 and the ranges of t' and x' are restricted to t' > 
and f(t',x') > 0. As we can check, it may be difficult to find such solution if we use 
comoving coordinates. Compare with its original form (|40|), which is really simple. 

The second family with q — 1 can be obtained by simply replacing the cosh (pt) 
functions by smh(pt) everywhere (so tanh(//t) — » coth(pt)). In this case the range for 
t is restricted to t > because there is an initial spacelike singularity at t — 0. 



5.3 n = 2 

Despite the fact that now there are two different subcases, we begin giving a common 
feature for them: we must put C3 = in order to have T linearly independent from 1 
and f 2 (see (0)). 

In the subcase (i) we take {l,T 2 ,Mg} as linearly independent functions. To this 
end, due to the definition of M 8 , we must also demand that c\ — C2/2 7^ 0. We impose 
then 

M 7 = c 71 T 2 + c 72 M 8 , 

which must be an identity, i.e. multiplying this expression by c\ — c 2 /2 we get a 
linear relation between the three linearly independent functions and thus the constant 
coefficients must vanish, providing the following constraints on the contants 

(n j- M c 2 + 2ci 

C71 = c 2 2ci + -c 2 , c 72 = — , a = c 72 J. 

V 2 J c 2 - 2ci 
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The rest of the qa constants are given by cli = {c 2 , c|, 0, cic 2 , 0, 0}, cl 2 = 0, cgi 



0,Cg2 = 1, while hi = 0. Regarding the t-functions, it only remains the relation (22), 
which gives the equation for T(t): 



Br 2, — 1C 

T 2 (2 Cl + c 2 ) + 2T(C 2 - 2[3c 72 ) + TT 2 c\ + 2T 2 (^ i - c 2 C 2 ) 

2ci - c 2 

4 2 , (/3 2 c 72 + /5C 2 + B) 



-2T*c Cl =4 

2ci - c 2 

The system for the x-functions (0), ( p)D and flT5[ ) is given by 

^3 + ^3^2 = 0) 

N 7 N 8 = B, 

6 

5> L1 iV L + c 71 iV 7 = d, 

L=l 

c 72 N 7 + N S =C 2 . 

This system contains four equations for three unknowns. Thus, taking into account 
that we must impose ^ to avoid a third isometry, the compatibilization of the 
system gives raise to a particular family of solutions, apart from another solutions 
which fall into the previous cases n = 0,1. After redefining some constants, the line 
element reads 

ds 2 = e 2 ^ +2ciT W (-dt 2 + dx 2 ) + e C2T « (e 2l/x dy 2 + e~ 2vx dz 2 ) , (41) 

where p and v are constants, c 2 7^ 0, c 2 + 2c\ 7^ and T(t) satisfies 

1 

2(c 2 T2ci) 



T = — — ^— — f 2eVA - Ac 72 u 2 - c\t 2 



where e 2 = 1 and A = [c 7 \T 2 + 2c 72 v 2 ) 2 + 4/i 2 c 2 c 72 T 2 . This family does not contain a 
barotropic equation of state except for cases with more isometries. The pressure, the 
energy density and the rest of the variables involved in the fluid quantities are given 
by 



p = e - 2 ^ +c ^ — c 71 f 2 (2Cl " ° 2 



Va), 



c 2 + 2ci \ 2c 2 

-2(nx+ Cl T(t)) C 2 f-rja _ r. „2 ~ C 1 



p + p = e -^cu m>_±^ CnT ^ _ 2c 72 S - 2-eVA , £ = c 2 pT, 
c 2 + 2ci V c 2 / 

w = - (c 71 T 2 + 2c 72 z/ 2 - eVA) , w x = (c 71 f 2 + 2c 72 v 2 + ev^) , 

Z ^ ' ^C 72 ^ ' 

from where it is straightforward to find the pefect-fluid quantities using the expressions 
given in the subsection [4.5|, and realize that the region defined by A = does not 
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represent any physical singularity but a focal zone for the fluid congruence. The non- 
zero components of the Weyl tensor read 

xPo + ^4 = 2e' 2h pu, ij - ij 4 = -e- 2/ V(c 2 - 2 Cl )T, ^2 = ^ 2h [(c 2 - 2c x )f - 4z/ 2 

where fi = px + ciT(t), so the Petrov type is I. 

Although the Ae region is defined in principle by two inequalities involving the 
function T(t) (i.e. p + p > 0, Wq > 0), it can be shown that after imposing the former, 
it suffices the evaluation of the latter on JF (i.e. in p + p = 0), giving then only a 
condition on the constants involved. This happens because E never vanishes at T . 

For the subcase (ii) we take {1,T 2 ,M 7 } to be three linearly independent functions 
and impose the relation 

M 8 = c 81 t 2 . 

Following the same procedure as in the previous subcase, we must have now C2+2ci 7^ 0, 
C2 — 2ci = 0, (3 = 0, and cgi = 0, thus in fact, M 8 = identically. Therefore, equation 
(p2]) has the form 

C{T 2 + C 2 M 7 + B = 0, 

from where it follows that C\ = C2 = B = 0. The constants cli and cl2 are the same 
as in the subcase (i) and also C71 = 0, C72 = 1. The equation (|TjJ) for A = 2 reads then 

8 
i=l 

thus w 1 = and therefore the possible solutions belonging to this subcase are always 
separable in comoving coordinates. 

To sum up this subsection, in n = 2 there is only a family of solutions. Its line 
element is given by (0) and has no barotropic equation of state. 



6 Concluding remarks 

The assumption of separability of the metric functions in non-comoving coordinates has 
been shown to be a valuable tool for the obtaining of inhomogeneous exact solutions. 
By exploiting the coordinate interchange symmetry explained in section 2, a purely 
mathematical classification is put forward and provides a systematic way of getting 
new diagonal G2 on S% solutions. This has been explicitly used in the present work to 
construct several new solutions by means of the analysis of the simplest case (m=l). It 
arises a relationship between the different cases of the classification and some physical 
properties, such as the existence of a barotropic equation of state and its explicit form. 
Furthermore, the use of non-comoving coordinates shows explicitly that the space-times 
obtained by the imposition of a particular Segre type of the energy-momentum tensor 
(for instance, perfect fluid) may be extendible in general and with a varying algebraic 
type through the extension. 
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A The Einstein tensor and the functions M$(t) and 
N,(x) 



Here we give the explicit expressions for the Einstein tensor for the metric given in 
fll2]) computed in the natural orthonormal co-basis (§), and the functions Mj(t) and 
Ni(x) (i = 1 . . . 8) in terms of the metric functions defined upon equations ( |I4"D and 
(|15|)- The non-zero components of Einstein tensor for (|l2|) in the frame (Q) are: 

1-2 -2 3 



<Soo 
^ii 

<S()1 

S22 



-2/1 



-2/1 



-2/1 



5; 



33 



-2/i 



1 I 2 + ^ I 2 _J 3 _A 2 ~ 4 A 2 _A 3 + A 1 A 2 

1 •■ 3 • 2 • 2 • • ' 

A 1 A 2 + ^ A 2 — A 3 _J 2 _ ^ J 2 _J 3 + J 1 J 2 



^1^2 + ^2 [X'i - J ~ 2 T 3 X 3 

A i + 2 2 4 2 _ A 2 A 3 - A 3 + A : 



-t - \f 2 - -t 2 + T 2 T 3 + T 3 - T 3 2 



//!//!/ / / // 
X 1 + -X 2 + -X 2 + X 2 X 3 + X 3 + X 



-fi - -f 2 - -T 2 2 
1 2 2 4 



T 2 T 3 - T 3 



To 



The equation (IT3) is now used to eliminate T 3 and X 3 in what follows. There is some 
freedom in defining the M and N functions, and we could have chosen the definitions 
(in particular for M^t) and Nl{x) (L = 1 . . .6)) keeping the x «-> i symmetry; but, 
as we deal with the t-functions, we define Mt(t) as simpler as possible. From we 
define Mz,(t) and Nl(x) as 



M, /, • M 2 '2 • M 3 



T 3 , M 4 = T^, M 5 = TiTg, M 6 = T 2 T 3 , (42) 



AT, 



yl2 
" A 2 5 



iv 2 = - - \x' 2 



2^2 > 



x 



-4A?, 



2A ^2 (-^l - 1^2) ) 



X 



4^, 



3 / \ (43) 

X 6 = 4X 3 X( - |X 2 . 



Similarly, from the right hand side of (^), and using the explicit expressions for the 
two factors, we can define M 7 (t), M 8 (t), N 7 (x) and N 8 (x) by means of 



Soo + S22 

Sll — 5^22 

and such that ([X5|) holds. Thus, 



-2/: 

-2/: 



^M^+iVsOr)) 
MMg^+iV^x)) 



t x t 2 -\f 2 -lf z -f, + a, M % =f{f 2 -\f 2 -\f 2 +T 1 + /3, 



(44) 



N 7 = X[X< 2 - \X» - 2X£ - X'l -f3, N 8 = X[X' 2 - \X» - \X' 2 2 + X'{ - a, (45) 
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where a and (3 are two arbitrary constants. With the help of these constants we can 
always set, for instance, M7(£)+const.— > M 7 (t), that is, we can absorve the constants 
added to M7, M§ into themselves. This feature is very useful to simplify the expressions 
and the equations for Ni(x). 

B Some useful results 

Here we give some lemmas concerning sets of linearly independent functions which 
will be useful for relating the integers m and n defined in the section |j. The first two 
lemmas relate a set of functions {fi} (i — 1 . . . r), from (an interval of) 1R to H, with its 
derivatives (denoted by a dot). The rest give some relations of the same kind between 
the sets {fi} and {fifj} (the combinations of products of two fuctions from {fi}). We 
will present an example of each type of proof to give an outline of the procedures 
involved. 

Lemma 1: Let {fi} be a set of r linearly independent C 1 functions. Then, there are 
at least r — 1 linearly independent functions among the set {fi}- 

Proof. Suppose that there were only r — 2 linearly independent functions among {fi}. 
This would mean that we have two independent relations between them 

r r 
i=l i=l 

where <2j and hi are constants (or functions depending on another independent variable). 
Integrating these relations we would get 

i=i 1=1 
where A and B cannot vanish. Therefore 

j^iBa, - Abjfi = 0, 
i=i 

and thus Bai — Abi = because {fi} are r linearly independent functions. But then, 
the two relations in (|B|) would be linearly dependent, which is a contradiction. 

Lemma 2: {fi} is a set of r linearly independent functions if and only if {fi, 1} is a 
set of r + 1 linearly independent functions. 

Lemma 3: Let f\ and f% be two linearly independent functions such that the inter- 
section of their supports is non empty. Then, the set {fifj}= {ff, fi, /1/2} contains 
three linearly independent functions. 
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Let us note that the condition on the supports avoids /i/ 2 = (in all the interval of 
definition of {/«})• 

Proof. Suppose, on the contrary, that there existed a relation af\ + 6/| + cf 1 f 2 = 0. 
Dividing this relation by f\ (in its support), we would get a polynomial for /2//1. This 
would imply that a, b, and c are all equal to zero, in order to avoid the proportionality 
of f\ and f 2 . 

Theorem 1 Let / 1; f 2 , and f 3 be three linearly independent functions such that the 
intersection of any two of their supports is non empty. Then, the set {fifj} contains 
at least five linearly independent functions. 
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